Inverse Monte Carlo procedure for conformation determination of macromolecules

A novel numerical method for determining the conformational structure of macromolecules is applied to idealized biomacromolecules in solution. The method computes effective inter-residue interaction potentials solely from the corresponding radial distribution functions, such as would be obtained from experimental data. The interaction potentials generate conformational ensembles that reproduce thermodynamic properties of the macromolecule (mean energy and heat capacity) in addition to the target radial distribution functions. As an evaluation of its utility in structure determination, we apply the method to a homopolymer and a heteropolymer model of a three-helix bundle protein [Zhou, Y.; Karplus, M. Proc Natl Acad Sci USA 1997, 94, 14429; Zhou, Y. et al. J Chem Phys 1997, 107, 10691] at various thermodynamic state points, including the ordered globule, disordered globule, and random coil states.